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Abstract 

Among the discrete evolution equations describing a quantum sys- 
tem Tis undergoing repeated quantum interactions with a chain of 
exterior systems, we study and characterize those which are directed 
by classical random variables in M. N . The characterization we obtain 
is entirely algebraical in terms of the unitary operator driving the el- 
ementary interaction. We show that the solutions of these equations 
are then random walks on the group U{7io) of unitary operators on 
Ho. 

1 Introduction 

In the article |APj Attal and Pautrat have explored the Hamiltonian descrip- 
tion of a quantum system undergoing repeated interactions with a chain of 
quantum systems. They have shown that these "deterministic" dynamics 
give rise to quantum stochastic differential equations in the continuous limit. 
Since that result, some interest has been found in the repeated quantum 
interaction model in itself (cf [A7T] . [Al2] . jBJMl] . [BJM2] , [BJM3] ) and 
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several physical works are in progress on that subject (for example jAKPj ) . 
These repeated interaction models are interesting for several reasons: 

- they provide a quantum dynamics which is at the same time Hamiltonian 
and Markovian, 

- they allow to implement easily the dissipation for a quantum system, in 
particular they are practical models for simulation. 

The probabilistic nature of the continuous limit found by Attal and Pautrat 
is not due to the passage to the limit, it is already built-in the Hamiltonian 
dynamics of repeated quantum interactions (it is actually built-in the axioms 
of quantum mechanics). 

The evolution equations decribing the repeated quantum interactions are 
purely deterministic but they already show up terms which can be interpreted 
as "discrete time quantum noises". The point with these discrete quantum 
noises is that sometimes they may give rise to classical noises. That is, some 
linear combinations of these quantum noises happen to be mutually commut- 
ing families of Hermitian operators, hence they simultaneously diagonalize 
and they can be represented as classical stochastic processes. 

In the other cases, that is, with different combinations of the quantum 
noises, no classical process emerges and the dynamics of repeated quantum 
interactions is purely quantum. 

The aim of the article is to explore the case when the dynamics is clas- 
sically driven. We characterise algebraically, on the Hamiltonian, the case 
when the dynamics is classical. 

The article is organised as follows. We first (Section EJ) present the phys- 
ical and mathematical setups for the repeated quantum interactions. In Sec- 
tion [3] we introduce the basic algebric tool: the obtuse random walk which are 
an appropriate "basis" of random walks adapted to this language. We then 
explore and characterise the unitary random walks which emerge classically 
from the repeated quantum interactions (Section HI). We finally specialize 
in Section [5] our result to the one dimensional case which alreay shows up a 
non-trivial structure. 

2 Repeated Quantum Interactions 
2.1 The Physical Model 

Repeated quantum interaction models are physical models developed by At- 
tal and Pautrat in [APJ which consist in describing the Hamiltonian dynamics 
of a quantum system undergoing a sequence of interactions with an environ- 
ment made of a chain of identical systems. These models were developed 



2 



for they furnish a toy model for a quantum dissipative system, they are at 
the same time Hamiltonian and Markovian, they spontaneously give rise to 
quantum stochastic differential equations in the continuous time limit. Let us 
describe precisely the physical and the mathematical setup of these models. 

We consider a reference quantum system with state space Ho, which we 
shall call the small system (even if it is not that small!). Another system 
He, called the environment is made up of a chain of identical copies of a 
quantum system H, that is, 

H E = (g)H 

nGN* 

where the countable tensor product is understood in a sense that we shall 
make precise later. 

The dynamics in between Ho and He is driven as follows. The small 
system Ho interacts with the first copy H of the chain during an interval 
[0, h] of time and following an Hamiltonian H on Ho <S> H. That is, the two 
systems evolve together following the unitary operator 

U = e~ ihH . 

After this first interaction, the small system Ho stops interacting with the 
first copy and starts an interaction with the second copy which was left un- 
changed until then. This second interaction follows the same unitary operator 
U . And so on, the small system Ho interacts repeatedly with the elements 
of the chain one after the other, following the same unitary evolution U . 

Let us give a mathematical setup to this repeated quantum interaction 
model. 

2.2 The Mathematical Setup 

Let Ho and H be two separable Hilbert spaces (in the following, for our prob- 
abilistic interpretations the space H will be choosen to be finite dimensional). 
We choose a fixed orthonormal basis {X n ; n G M U {0}} where M = N* or 
{1, . . . , N} depending on wether H is infinite dimensional or not (note the 
particular role played by the vector X° in our notation). We consider the 
Hilbert space 

where this countable tensor product is understood with respect to the stabi- 
lizing sequence (X°) nGN ». This is to say that an orthonormal basis of T<E> is 
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made of the vectors 

X a = <g) xt 

neN* 

where a = (i n )neN* runs over the set V of all sequences in A/" PI {0} with only 
a finite number of terms different of 0. 

Let U be a fixed unitary operator on H$®H. We denote by U n the natural 
ampliation of U to Ho <8> where C/ n acts as [/ on the tensor product of 
Ho and the n-th copy of H and U acts as the identity of the other copies of 
H. In our physical model, the operator U n is the unitary operator expressing 
the result of the n-th interaction. We also define 

v n = u n u n _ 1 ...u 1 , 

with the convention Vq = I. Physically, V n is clearly the unitary operator 
expressing the transformation of the whole system after the n first interac- 
tions. 

Define the elementary operators a*, i,j G Af H {0} on H by 

a)X k = 6, X X > . 

We denote by a*-(n) their natural ampliation to T$ acting on the n-th copy 
of H. only. That is, if a — (? n ) Tte N* 

a)(n)X a = 5i,i(n)X a \ {in}u{jn} . 

One can easily prove (in the finite dimensional case this is obvious, in the 
infinite dimensional case it is an exercise) that U can always be written as 

i,jeATu{0} 

for some bounded operators C/j on Ho such that: 
- the series above is strongly convergent, 

With this representation for U, it is clear that the operator U n , repre- 
senting the n-th interaction, is given by 

U n = Uj®a)(n). 
j,jeA/u{0} 
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With these notations, the sequence (V n ) of unitary operators describing the 
n first repeated interactions can be represented as follows: 

Ki+i = u n+ i v n 

= U)®a){n + l)V n . 

But, inductively, the operator V n acts only on the n first sites of the chain 
T$, whereas the operators a*(n + 1) acts on the (n + l)-th site only hence 
they commute. In the following, we shall drop the ® symbols, identifying 
operators like a* (n + 1) with I no ® a* (n + 1). This gives finally 

V n+1 = U]V n a){n + l). (1) 

i,i &Vu{0} 

In Quantum Probability Theory, the operators a* (n) have a particular 
interpretation, they are discrete-time quantum noises, they describe the dif- 
ferent types of basic innovations than can be brought by the environment 
when interacting with the small system. See |Atj for complete details on 
that theory, the understanding of which is not necessary here. 

The only important point to understand at that stage is the following. 
In some cases the above equation ([T]) corresponds to an equation divren by a 
classical noise, i.e. driven by a random walk. This is what we shall describe 
in the next section. 

3 Classical Random Walks 

In order to understand the link that may exist between the discrete-time 
quantum noises a* and classical random walk, one needs to pass through 
a particular family of random walks, the obtuse random walks. Defined by 
Attal and Emery in |A-Ej . these random walks constitute a kind of basis of 
all the random walks in M . Let us describe them. 

3.1 Obtuse Random Walks in 

Let X be a random variable in taking A/" + 1 values Vq, vn with respec- 
tive probabilities p , ...,Pn such that p« > 0, Wi G {0,1, N}. The canonical 
space of X is the triple (A, A, P), where A = {0, 1, A^}, A is the a— field 
of subsets of A and P is the probability measure given by P({i}) = Pi- Hence 
for all i G {0, 1, A^} we have X(i) = v { and P(X = v { ) = P({i}) = Pi- 
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We say that such a random variable X is centered if E(X) = (as a vector 
of K^). We say that X is normalized if Cot> (X) = / (as a iV x iV- matrix). 

Let us denote by X 1 , ...,X the coordinates of X in the canonical basis 
of R and define the random variable X° on (A, .4., P) given by X°(i) = 
1, Vz G A. Let us introduce the random variables X % defined by 

X i (j) = y/PjX^j), 

for all i, j G {0, 1, N}. We then have the following easy characterization 

(cfra. 

Proposition 3.1 The following assertions are equivalent: 

1 ) The random variable X is centered and normalized, 

2) The family vq, ...,Vn of values of X satisfies < Vi , Vj >= —1, for all 
i 7^ j and the probabilities pi 's are given by 

Pi = : — n n 2 , for alii G {0,1,..., N} , 



3) The matrix (X°, X 1 , X ) is unitary. 
A family of N + 1 vectors in ~R N satisfying the above condition 

< Vi, Vj >= -1 , 

for all % ^ j, is called an obtuse system in [A-E] . Hence, a random variable 
X satisfying one of the above condition is called an obtuse random variable. 

Note that, as a corollary of the above proposition, the random variables 
X°, X 1 , ...,X N are linearly independent and hence they form an orthonormal 
basis of L 2 (A, A, P). In particular, for every i,j G {1, . . . , N} the random 
variable X l X^ can be decomposed into 

N 

x { x j = T k xk 

for some real coefficients Tjf . The familly of such coefficients forms a so- 
called 3-tensor, that is they are the coordinates of a linear mapping T from 
R N to Mjv(R). 

We say that a 3-tensor T is ses qui- symmetric if the two following assump- 
tions are satisfied: 
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i) (i, j, k) i — > Tl J is symmetric, 

ii) (i, j, I, m) i — > J2k=i T k T k m + SiAm is symmetric. 

Using the commutativity and the associativity of the product X l X^ it is 
easy to prove the following (cf [A-E] ) . 

Theorem 3.2 If X is a centered and normalized random variable in M , 
taking exactly N + 1 values, then there exists a ses qui- symmetric 3-tensor T 
such that 



In the following, by an obtuse random walk we mean a sequence (X p ) peN 
of independent copies of a given obtuse random variable X. Actually, the 
random walk is the sequence made of the partial sums Yl p < n but we shall 
not make any disctinctions between the two processes in the terminology. 

3.2 More General Random Variables 

We claimed above that obtuse random variables are a kind of basis for the 
random variables in R^ in general. Let us make precise here what we mean 
by that. 

First of all, a remark on the number AT + 1 of values attached to X in 
Mr . If one had asked that X takes less than AT + 1 values in R^ (k, say)and 
be centered and normalized too, it is not difficult to show that X is actually 
taking values on a proper subspace of M N , with dimension k — 1. For example, 
a centered, normalized random variable in R 2 which takes only two different 
values, is living on a line. 

Now, if Y is a random variable in R^ taking k different possible values 
Wi, . . . , Wk, with probability pi, . . . ,p k and k > n+1. Consider an obtuse ran- 
dom variable X in R fc_1 taking values v±, . . . , with the same probabilities 
Pi, . . . ,Pk as those of Y. We have seen that the coordinate random variables 
X 1 , . . . ,X k ~ 1 , together with the deterministic random variable X°, form an 
orthonormal basis of L 2 (A, A, P) we can represent each of the coordinates of 



j=0 

Hence we have a simple representation of Y in terms of a given obtuse random 
variable X. 



X®X = I + T(X). 



Y as 



fe-i 
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3.3 Connecting With the Discrete Quantum Noises 

The obtuse random walks admit a very simple and natural representation in 
terms of the operators a* (n) defined in Section 12.21 

Let X be an obtuse random variable in M. N . On the product space 
(A N , A® N , P® N ) we define a sequence (X p ) pgN of independent, identically 
distributed, random variables, each with the same law as X. 

Consider the space T$(X) = L 2 (A N , A m , P m ) and the random vari- 
ables 

Xa= J] 

where A is any sequence in {0, 1, iV} with only finitely many terms differ- 
ent from 0. 

The following result is also easy to prove (cf [At] ) . 

Proposition 3.3 The random variables Xa, where A runs over the sequences 
in {0, 1, ...,N} with only finitely many terms different from 0, form an or- 
thonormal basis o/T$(X). 

In particular we see that there exists a very natural Hilbert space isomor- 
phism between the space T$(X) and the chain T$ constructed in Section 
12.21 over the space 7i = C N+1 . Regarding this isosmorphism, one can con- 
sider the operator M.x^p) of multiplication by the random variable X p on 
T$(X). This self-adjoint operator contains all the probabilistic information 
associated to the random variable X p , it admits the same functional calculus, 
etc ... it is the actual representant of the random variable X l p in this Hilbert 
space setup. 

As each of the probabilistic space T$(X) are made isomorphic to T$ we 
can naturally wonder what happens to the operators M.xHp) through this 
identification. The answer is surprisingly simple (cf [At] ) . 

Theorem 3.4 Let X be an obtuse random variable in M. N and let (X p ) pe N 
be the associated random walk on the canonical space T$(X). Let T be the 
ses qui- symmetric 3-tensor associated to X. If we denote by U the natural 
unitary isomorphism from T$(X) to T$ , then for all p G N, i G {1, N} 
we have 

N 

UM x dJ* = a°(p) + 4(p) + ]T T? l al(p) . 

3,1=1 

Here we are! By a simple linear combination of the basic matrices a*(p) 
one can reproduce any random variable on ~R N . 
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Coming back to the evolution equation (JT]), we see basically two different 
cases may appear. 

First case: the coefficients C/j of the basic unitary matrix U are such that 
the equation CD reduces to something like 

N 

V n+1 = AV n + J2 BiV n Mxi ■ 

i=l 

This means that this operator-valued evolution equation, when transported 
back to T$(X) is an operator-valued (actually unitary operator-valued) 
equation driven by a random walk (X p ) gN . It is a random walk on U(N). 

Second case: there is no such arrangement in the equation (TJ this means it 
is purely quantum, it cannot be expressed via classical noises, only quantum 
noises. 

Our aim, in the rest of the article is to characterize completely those 
unitary operators U which give rise to a classically driven evolution (first 
case) . 



4 Random Walks on U(Hq 



In this section we work on the state space 

T$ = (g) C N+1 . 

neN* 

We consider a fixed obtuse random variable X, with values v\, . . . , v n and 
with associed 3-tensor T. We identify the operator 

N 

«°(P) + 4(P) + E^V(P) 

3,1=1 

with the random variable X % v and we denote it by X*, instead of M.xy Recall 
that X® is the constant random variable equal to 1, hence as a multiplication 
operator on T$ it coincides with the identity operator I. 

In the following we extend the coefficients of the 3-tensor T to the set 
{0, 1, . . . , N}. This extension is achieved by assigning the following values: 

rpij rpiO rpOi c 

1 - lj - lj - Oij . 

With that extension, the second sesqui-symmetric relation for T is written 
simply 



N 

n) (i, j, I, m) i — > T^T 1 ™ is symmetric. 

fc=0 
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Recall the discrete time evolution equation §T§ associated to the repeated 
quantum interactions: 

N 

V n+1 = U]V n a){n + l), 

i,j=0 

with the convention Vq = I. 

Proposition 4.1 The discrete-time evolution equation |7]) can be written as 

N 

Vn+1 = £ VnX % n+ l , 
i=0 

for some operators B k on Ho, if and only if the coefficients t/j are of the 
form 

N 

0J = £ ■ (2) 

k=0 

Proof Let us prove first the sufficient direction. If U is of the form (|2j) then 

N 

V n+1 = + 

i,j=0 

N N 

= U°V n a (n + l) + YUoV n a (n + l) + J2u?V n a i (n + l) + 

i=l i=l 

N 

+ Y U]V n a)(n + l). 
The relation (J2]) implies in particular Uq = B and Uf — Uq — Bi. This gives 

N 

V n+1 = B V n a (n+l) + ^B i VM(n+l) + a° i (n + l)) + 

i=i 

N N N 

+ £ J2 Tl j B k V n a){n + 1) + £ B V n aj(n + 1) 

fc=l i,j=l i=0 

N N 

= B V n + B kVn [4(n + 1) + a°(n + 1) + £ 7fa}(n + 1)] 

k=l i,j=l 
N 

= B V n + ^ -BfcKi -^"n+l 

fc=l 

N 
k=0 
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This gives the requiered result in one direction. The converse is easy to prove 
by reversing all the arguments above. □ 

Now, consider the operators 

N 
i=0 

with the convention v® = 1, for all k G {0,1, N}. Our purpose in the 
sequel is to prove that these operators are unitary if and only if the evolution 
operator U is unitary. Here is the first step. 

Proposition 4.2 If U a unitary operator, then for all I G {0, 1,...,N} the 
operator W\ is unitary. 

Proof We have 

N 

Wiwt = E M B &) ■ 

i,j=0 

But expressing coordinate-wise the relation §2§, we have 

N 

M = E w ■ 

m=0 

Hence, we get 

N 

ij, m=0 

N N 

= E ^ n EW fl i 

j,m=0 i=0 
N 

= E V TU 3 JJT 

j,m=0 
N N 

= E^E^f) 

m=0 j=0 
N N 

= Z> m (£*«* J ) 

m=0 j=0 
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This completes the proof. □ 

Now, our aim is to prove the converse of Proposition 14.21 In order to 
achieve this, we need to express the coefficients Uj of U in terms of the 
operators Wj's. This is the aim of the following two lemmas. 

Lemma 4.3 For all i G {0, 1, N} we have 

N 



1=0 

Proof We have 

N N N 

E?" 1 '?^ = E^(E^) 

1=0 1=0 j=0 

N N 

= E^-(EmV) 

3=0 1=0 
N 

= ^2 B j E(X i X j ) 

j=0 
N 

= E B 3 Si 3 = Bi ■ 
3=0 

This ends the proof. □ 
Lemma 4.4 For all I, k G {0, 1, N} we have 

A' 



8=0 

Proof Recall that we have 



N 



uf = E T f B i 

3=0 



and 



N 



( 3 ) 

3=0 
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By using Lemma [4.31 and relation ([3]), we get 

N 

U? = Y. I';'' j ' <■■"'> 

i,j=0 
N N 

i=0 j=0 
N 

= ^ PiV$v\Wi. 

□ 

As a corollary of the two above lemmas, we prove the following. 

Proposition 4.5 If all the operators Wi, fori e {0,1, N} , are unitary 
then the operator U is unitary. 

Proof We have 

TV 

£(tfj)(tf£)* = E PiPj^vlvfWiW* 

k=0 i,j,k=0 

N N 

= p* (v^vr i + e ''■i>J r ' i 1 r '>i' 

i,k=0 i,j,k=0,ij^j 
N 

i=0 

N N 

+ E PiPj(J2 v i v 1) v i v 7 w i w ! 

i,j=0,iy^j k=0 

N 
i=0 

N 

+ E ^;(<^%>+i)^ m ^/- 

But recall that, by Proposition ^. 11 we have Pi(||i>i| | 2 + 1) = 1 and < V{, Vj >= 
— 1 for all i j. Therefore we get 

N 

Y(U l Mr=E(X l X™)I = 5 ml I. 

k=0 
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We have prove the unitary character of U*. □ 
Alltogether we have proved the following result, which resumes all the 
results obtained above. 

Theorem 4.6 Let X be an obtuse random walk in M. N , with values v , . . . , vn, 
with probabilities Pq, . . . ,Pn and with 3-tensor T. Let (X p ) peN be its associ- 
ated obtuse random walk. Then the repeated quantum interaction evolution 
equation 

N 

K+i= ^t/jVU>+ 1 ) 

i,j=0 

takes the form 



N 

if and only if there exists unitary operators Wi, % G {0, . . . , iV} 7 on 7{q such 
that the coefficients C/j of U are of the form 



N 



uf = j2p* v i v i w *- 

In that case, the coefficients above are given by 

N 
1=0 

When the conditions above are satisfied, the evolution equation 

N 
k=0 



is, when seen in the space T<b(X), an operator- valued evolution equation, 
driven by a random walk. It is natural to wonder what king of stochastic 
process it gives rise to. 

Theorem 4.7 As a random sequence in U(7io), the solution of the equation 

N 

Vn+1 — ^2 BkVnXn+i 
k=0 



is an homogeneous Markov chain on U(N) (actually a standard random 
walk), described as follows: V — I almost surely and V n+ \ takes one of the 
values WiV n , % G {0, 1, N}, with respective probability pi, independently of 
V„. 
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Proof Assume V n is given, depending on the random variables Xi, . . . , X n 
only. Then the random variable X n+1 is independent and X l n+1 = vj, with 
probability p\. Therefore, with probability p\ we get 



V n+1 = ^B i vjV n = W l V n 



i=0 

This proves the result. □ 



5 The Case N = 1 

In order to illustrate the results of the previous section, we detail here the 
situation in the case N = 1. 

Consider the set Q = {0, 1} N , equiped with the cr-field T generated by 
finite cylinders. We denote by v n the coordinate mappings, for all n e N, 
that is z/ n (c<j) = uo{n). 

For p g]0, 1[ and q = 1 — p, we define the probability measure /i p on 

JF) which makes (z/„) ng N to be a sequence of independent, identically 
distributed, Bernoulli random variables with law pSi + qS . We denote by E p 
the expectation with respect to /i p . 

Define the random variables 



They satisfy E p [X n ] = and E P [X^] = 1, hence they are obtuse random 
variables in E. They take the two values t>o = \/ llv an d v i = ~^/p/l with 
respective probabilities p and q. 

The 3-tensor T associated to X is easy to determine. Indeed, one can 
easily check the following multiplication formula. 

Proposition 5.1 We have 

= 1 + C p X n , 

where c v = 

This means that the 3-tensor in this context, which is a constant, is T = c p . 

In this context also, note that the space T$(X) is the space L 2 (Q, J 7 , n P ), 
whereas the space T$ is ^C 2 . As an application of Theorem 13.41 the 
operator of multiplication by X n on T$(A) is represented on T$ as 

M x n = a°(n) + al(n) + c p a{(n) . 



s/pq 
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Here we are, we have put all the corresponding notations. We can apply 
Theorem 14.61 to this particular case. 

Theorem 5.2 Consider the obtuse random walk (X n ) neN on M., as described 
above. Then the repeated quantum interaction evolution equation 

N 

V n+1 = J2 u Fn4(n + l) 

i,j=0 

takes the form 

Ki+i = BqV„, + BiV n X n+ i 
if and only if there exist 2 unitary operators W and W\ on Hq such that 

JJ= ( pWo + qW, y/pq{W -Wi) 
~ Ky/miWo-Wx) qW + P W 1 

In that case, the coefficients Bi above are given by 

B = U° , B x = Ul = U l . 

The random sequence (V n ) n< - N is defined by V = I and 



V n 



+ 1 



Wo V n with probability p 
Wi V n with probability q . 
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